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Abstract 

This paper studies a prototype of inverse obstacle scattering problems whose gov- 
erning equation is the Helmholtz equation in two dimensions. An explicit method 
to extract information about the location and shape of unknown obstacles from 
the far field operator with a fixed wave number is given. The method is based 
on: an explicit construction of a modification of Mittag-LefHer's function via the 
Vekua transform and the study of the asymptotic behaviour; an explicit density in 
the Herglotz wave function that approximates the modification of Mittag-LefHer's 
function in the bounded domain surrounding unknown obstacles; a system of in- 
equalities derived from Kirsch's factorization formula of the far field operator. Then 
an indicator function which can be calculated from the far field operator acting on 
the density is introduced. It is shown that the asymptotic behaviour of the indicator 
function yields information about the visible part of the exterior of the obstacles. 
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1 Introduction 

This paper is concerned with developing an explicit analytical method for so-called inverse 
obstacle scattering problems at a fixed wave number. For the purpose we consider an 
inverse obstacle scattering problem in two dimensions in which the governing equation 
is given by the Helmholtz equation. The problem is to extract information about the 
location and shape of unknown sound-hard obstacles D embedded in a medium with 
constant acoustic speed and density, from the leading term of the asymptotic expansion 
of the reflected wave w at infinity which is caused by an incident plane wave e lkx ' d for 
infinitely many incident directions d E S 1 and a fixed wave number k > 0. This is a 
prototype of several inverse obstacle scattering problems of acoustic wave. 

More precisely we assume that: D C R 2 is open and R 2 \ D is connected; 3D is 
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Lipschitz. The reflected wave w is the unique solution of the scattering problem: 

Aw + k 2 w = in R 2 \ D, 

dw d ikx-d an 

— — = - — e lKxa ondD, 
av av 

lim Jr- l^r- — ikw] = 

r — >oo y (j r j 

where v is the unit outward normal relative to dD and r = \x\. This last condition is 
called the Sommerfeld radiation condition. 

It is well known that given ip <E S 1 the value w(rip) as r — > oo has the following form: 

e ikr / 1 \ 

W M = d,k) + o j . 

The coefficient Fr>((p; d, k) is called the far field pattern of w. 
The operator F : L 2 (S l ) — >■ L 2 (S l ) given by the formula 

F D g{ V )= f F D (<p;d,k)g(d)da(d), g e L 2 {S l ), 

is called the far field operator. It is welll known that the far field operator for a fixed 
k uniquelly dertmines the obstacles [19]. In this paper we consider: how to extract 
information about the location and shape of D from the far field operator or its partial 
knowledge at fixed k. 

In [9, 10] the author established the reconstruction formula of D itself from the far field 
operator. The formula consists of two parts: a relationship between a suitable Dirichlet- 
to- Neumann map on the boundary of a domain that contains D and the far field operator; 
application of the probe method introduced by the author [8] to the Dirichlet-to-Neumann 
map. 

In [24] Potthast gave a reconstruction procedure that he calls the singular sources 
method. The method yielded a way from the far field operator to a scattered field outside 
unknown obstacles which was exerted by a point source outside the obstacles and blows 
up on the boundary of the obstacles. 

Kirsch gave two types of reconstruction formulae of D in [17, 18]. The idea behind the 
formulae is called the factorization method since the formulae are based on a factorization 
formula of the far field operator. In particular, in [18] he made use of the quadratic form 

(F D g,g)L2(si) = J sl (F D g)(^)g(<p)do-((p) 

acting on densities g G L 2 (S' 1 ) to introduce his indicator function. It is defined by 

K(x) = M{\(F D g,g)\\ge L 2 (S r ), ($ x , g)^) = 1} 
where & x {(p) = e~ tkx ' v ', ip G S 1 . He established the one line formula 

D = {x G R 2 | K(x) > 0}. 
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For applications of his method to obstacles with other boundary conditions see [5, 6]. 

In [11] in three dimensions the author gave an extraction formula of the convex hull 
of D with a constraint on the Gaussian curvature of dD from a Dirichlet-to-Neumann 
map calculated from the far field operator. See also [12] for the sound-soft obstacles. It 
is an application of the enclosure method introduced by the author [12] and based on the 
asymptotic behaviour of the function 

having large parameter r where both $ and ft 1 - are unit vectors and perpendicular to each 
other. This function satisfies the Helmholtz equation Av + k 2 v = in the whole space 
and divides the whole space into two parts: if x ■ $ > t, then e~ Tt \v | — > oo as r — > oo; 
if x • $ < t, then e~ Tt \v | — ► as r — > oo. The indicator function introduced in [11] tells 
us whether given t the half space x ■ d > t touches unknown obstacles. 

The aim of this paper is to generalize this result by introducing another indicator 
function which is given by the form (F D g, g) acting on explicit densities g on S* 1 and tells 
us whether a given cone touches unknown obstacles. 



1.1 Statement of the main result and a corollary 

In this paper we identify the point •d = (i?i,i?2) £ •S' 1 with the complex number + ii? 2 
and denote it by the same symbol i?. 

Definition 1.1. Given n > 1, N — 1, • • •, s > and (y, uj) G Br x S l define the indicator 
function 



IN 



L 2 (5!) 



where 



i/n ( t \ ^ T(m+1) (suj^Y x 

^ r bd (i.i) 

Let < a < 1. Let C y (u, ira/2) denote the interior of the cone about u of opening angle 
7tq;/2 with vertex at y: 

C y (uj, ira/2) = {x e R 2 | (x — y) ■ uo > \x — y\ cos (net/ 2)}. 

The following theorem is the main result of this paper. 

Theorem 1.1. Let k 2 be not a Neumann eigenvalue of —A in D. Assume that dD is 
C 2 and that D is contained in the open disc Br radius R centered at the origin. Let 70 
be the unique positive solution of the equation logt + t/e = 0. Let 7 satisfy < 7 < 70. 
Let n > 1. Let {s(N)} N=1 ... be an arbitrary sequence of positive numbers satisfying, as 
N — > 00 

(Rs(N)) n = In + 0(1). 

e 

Then, given (y, u) e Br x S* 1 we have: 

ifC y (LU,n/2n)nD = 0, then lim^oc \I^M n ))n\ = 0; 
ifC y (u,n/2n) nD^I, then lim^oo \ll£ u) (s(N)) N \ = 00. 
Theorem 1.1 is a direct consequence of two lemmas below. 
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Lemma 1.1. There exists a positive constant C such that, for all g e L 2 (S' 1 ) 

\\Hg\dD\\ 2 H i/2(dD) < \( F Dg,g)msi)\ < C\\Hg\ dD \\ 2 Hl , 2{dD) 

where 

Hg{x)= e ikx -v g (ip)do-(tp), xEB 2 
J s 

and is called the Herglotz wave function with the density g. 

Lemma 1.2. Given (y,ou) G B R x S* 1 we have: 

if C y (uj,7r/2n) n D = 0, t/ien lim^-^oc | )( • ! s(iV), k,u) N \\ H i( D ) = 0; 
ifC y (u,7r/2n) n/)^0, toen lim^^^ H-ZfyJ^, } ( • ; s(N),k,uj) N \ dD \\ L 2 (dD) = oo. 

Lemma 1.1 has been pointed out in [5]. It is a corollary of a factorization formula in 
[17] and a coerciveness of an operator in the formula. A known proof of the coerciveness 
is given by a contradiction argument (cf. Lemma 4.2 in [6]) and therefore not direct. It 
seems that at the present time, there is no direct proof of this fact. 

Lemma 1.2 follows from Corollary 2.1 and (3.16) in Sections 2 and 3, respectively. 

So from Theorem 1.1 what information about unknown obstacles was extracted? To 
answer precisely we formulate the visible part of B R \ D. 

Definition 1.2. We say that a point y in B R \D is visible if the point y can be connected 
with infinity by a straight line that started at y and goes to infinity without intersecting 
D. We denote by V(B R \ D) the set of all points in B R \ D that are visible. We call this 
set the visible part of B R \ D. 

It is easy to see that the point y in B R \ D belongs to the visible part of B R \ D if 
and only if there exist n and u E S 1 such that {C y (u, ira/2)} fl D = with a = 1/n. 
The set V(B R \ D) is a non empty open set of B R \ D. If D is convex, then we have 
V(B R \ D) — B R \ D by a separation theorem. Needless to say, in general this is not 
true, however, the complement of the visible part of B R \ D gives an estimation of D from 
above. 

The next theorem tells us that the asymptotic behaviour of the indicator function 
J^^(s(A r )) A r as iV — > oo for all n and (y, uo) e B R x S* 1 uniquely determines the visible 
part of B R \ D except for a thin set. 

Corollary 1.1. Let D\ and D 2 be two obstacles such that: k 2 be not a Neumann eigen- 
value of —A in Dj; dDj is C 2 and that Dj is contained in the open disc B R radius R 
centered at the origin. Assume that, for each fixed n and (y, u) e B R x S* 1 we have 

Then V(B R \ D x ) \ dD 2 = V(B R \ D 2 ) \ dD,. 

This is derived from Theorem 1.1 as follows. It suffices to prove V(B R \ Di) \ dD 2 C 
V{B R \TJ 2 ). L et ye V{B R \~D 1 )\ dD 2 . Then there exist n and u e S 1 such that 
{C y (uj, na/2)} n D 1 = with a = 1/n. By Theorem 1.1 we have 
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Since y does not belong to dD 2 , it belongs to D 2 or B R \ D 2 . If y E D 2 , then from 
Theorem 1.1 one has 

A combination of this and assumption gives 

This is a contradiction. So y has to be in the set Br \ D 2 . If {C v (oj, na/2)} C\ D 2 ®, 
then, from Theorem 1.1 we again obtain (1.2) and the same contradiction as above. Thus 
{Cyiyj, na/2)} fl D 2 = 0. Therefore if one chooses a larger n' than n, then one gets 
{Cylco, ira'/2)} nD 2 = ® with a' = 1/n'. This means that y E V{B R \ D 2 ). 



1.2 A brief explanation of the idea 



Here we give a brief explanation of the origin of the density g}y n w \{ • ', s, fc)jv- 

Finding the density is closely related to Mittag-Leffler's function E a (z) which is an 
entire function and defined by the formula 

oo . m 

E a {z) = E 



m=0 



T(am + 1) 



where a is a parameter and satisfies < a < 1. The function E a (r(x\ + ix 2 )) of in- 
dependent variables x = (xi,x 2 ) with parameter < r < oo is harmonic in the whole 
plane. This function divides the whole plane into two parts as r — > oo: in a sector 
it is exponentially growing; outside the sector decaying algebraically. In [10] we applied 
this property of the harmonic function to an inverse boundary value problem for an el- 
liptic equation V • 7VM = with a discontinues coefficient 7 which is a special, however, 
very important version of the Calderon problem [3] and a continuum model of electrical 
impedance tomography. 

In Section 2 we modify this harmonic function by using the Vekua transform [26, 
27] (see also [2, 4]) which transforms given solution of the Laplace equation in R 2 into 
that of the Helmholtz equation Av + k 2 v = in R 2 . Using the solution obtained by the 
transform, we define a special solution with a large parameter s > of the Helmholtz 
equation which is denoted by E a (x;s,k,cu). In particular, the function Ei(x;s,k,u) is 
the Vekua transform of the harmonic function 



exp 1 1 (wi - iu 2 ) (xi + ix 2 ) j 



where uj = (cui,u! 2 ) E S 1 . We show that the function E a (x; s, k,u) has the asymptotic 
behaviour as s — > 00 similar to that of the original Mittag-Leffler's function: 

if x E Co (a;, na/2), then lim s >OQ \E a (x; s, k,u)\ = 00; 

if x E R 2 \ Cq(uj, ira/2), then lim s \E a (x; s, k,u)\ = 0. 

In Section 3 we establish the relationship between the density g^^i • s, Jc)n and the 
function E a (x — y; s, k, uj) for a — 1/n and y E Br. 

#<4 A L)( • 5 s(N), k) N (x) w E 1/n (x - y; s(N), k,u), x EB R 
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as N — > oo. Thus one can say that g]^^ ( • ; s, k)^ and s(N) are chosen in such a way that 
the corresponding Herglotz wave function approximates a modification of Mittag-Leffler's 
function. 

It should be pointed out that the result in [15] is closely related to the construction 
of the density. Therein the author considered the case when D is polygonal. This means 
that D has the expression D = Di Li ■ ■ ■ Li D rn with 1 < m < oo where D 1: ■ ■ ■ , D m are 
simply connected open sets, polygons and satisfy Dj fl Dj> = for j ^ f. 
The observation data are given by Fd{ • ;d,k) for fixed d and k provided we know the disc 
Br that contains D. Using the enclosure method [13], we established a direct extraction 
formula of the convex hull of D from the quantity 

i F D (-<p;d, k)g(tp)da(tp) 

for some explicit densities g independent of D. The one of key points is the choice of the 
densities. Those are chosen in such a way that 

Hg(x) « e *-(™W^Wu,^ x e s R 

where u 1 - = (— u^^i)- However, to get more than convex hull of unknown obstacles the 
function in the right hand side is not enough. In this paper we give explicitly the desired 
function by using the idea of the Vekua transform and Mittag-Leffler's function. 

Finally we point out that there are other approaches with a single incident plane wave: 
the point source method [23], the no response test [21], the range test [25] and the notion 
of the scattering support [20]. 



2 Modified Mittag-Leffler's function 

In this section we introduce a modification of E a (r(xi + ix 2 )) with < r < oo that 
satisfies the Helmholtz equation Au + k 2 u = in R 2 and study its asymptotic behaviour 
as r — > oo. 

The Bessel function of order m — 0, 1, • • • is given by the formula 



j (t) = r-Yy ( y ('- 



% (m + n)\n\ \2 
Definition 2.1. Let k > and < a < 1. Define 

Using the well known inequality (see Ex. 9.6, p.59 of [22]) 

ft \ m 1 

\Jm(t)\<{^) ^tGR, (2.1) 

one knows that E*(x; r) is well defined and satisfies \E*(x; t)\ < E a {r\x\). 

The idea behind Definition 2.1 is the following. Let x = (rcos#, rsin#). One has 



oo T 

E a (r(x 1 + ix 2 )) = ]T — — 



m 

m Amd . 



q r(am + 1) 



r e 
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oo m, r / 9 \ m 1 

£ ° ( ^ ) = ro ? r(^TT){(l) 

Therefore, by replacing r m in the expansion of E a (r(xi + ix 2 )) with (2/k) m m\J rn (kr) one 
obtains E*(x; r). From (2.1) one knows that the absolute value of (2/k) m m\J m (kr) is not 
greater than r m and 

/2\ m 

f -J m\J m (kr) ~ r m 

as r — > 0. In particular, we have E^(0;t) = 1. 

Using the change of variables w = y/l — t and the formula 

k\x\ [\l - w 2 ) m JAk\x\w)dw = 1 - ( m\J m (k\x\), 
Jo \k\x\J 

one can see that the function E^(x; r) has the integral representation given by the formula 

k I x I f i dt 
E k a { X] t) = E a (r( Xl + ix 2 )) - -LI / £ Q (<rt(m + ix 2 ))J 1 {k\x\VT^t)^=. (2.2) 

The integral transform 

— v(x) ; — - / V (tx) Ji(k\x\y/l — t) = 

2 Jo v 1 — t 

is called the Vekua transform of the function v(x) into the function u(x). This transforms 
given solution of the Laplace equation in R 2 into that of the Helmholtz equation Au + 
k 2 u = in R 2 . The formula (2.2) says that E*(x; r) is the Vekua transform of E a (r(x\ + 
ix 2 )) and therefore satisfies the Helmholtz equation. 

In this section we show that E^(x;r) as r — > oo has the almost same asymptotic be- 
haviour as E a (r(xi + ix 2 )). 
In this paper, for convenience we introduce 



/2\ m 

Jm(t) = (jj m\J m (t) 



From (2.1) we have \J m (t)\ < 1. In what follows we make use of this inequality frequently. 

Let f(z) be an arbitrary entire function of independent variable z = X\ + ix 2 . Let 
u(x; r) denote the Vekua transform of f(r(xi + ix 2 )). u(x; r) takes the form 

u{x-t) = f(r(x 1 + ix 2 )) - f-j^J J f(rt(x 1 +tx 2 ))J 1 (k\x\VT :r t)dt. (2.3) 

In this section we write C{ira/2) = C ((l, 0), ira/2). 

The following is useful for the treatment of E^(x;t) outside the cone C(na/2). See 
Appendix for the proof. 

Lemma 2.1. One can write u(x; r) anci the partial derivatives as: 

( 2 \ 2 1 /" r 

Hx;r) -/CKxi + iara))} = -Ji(fc|a;|)- / f(w(x 1 + ix 2 ))dw + R(x; r) (2.4) 

r Jo 
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where R(x; r) satisfies 



\R{x-r)\<U^f \ ^£w\f(w( Xl + ix 2 ))\dw; (2.5) 



2 \ f d , N ., , Pl , , \\\ V~ X Ji 



x 



u(x; t) - tV /'(r(xi + ix 2 )) ) = — f(r(xi + ix 2 )) 



y /.• x y [ r'Ar, ' J xi + ix 2 

+ JH^^MM + (^) 2 Xj J 2 (/c|x|)| - [ T f(w( Xl + tx 2 ))dw + Rjix; r) 



(2.6) 



xi — ix 2 
where Rj(x; r) satisfies 

|^(x;r)|<- f-Lij x 



^ jfV|/>(xi+ix 2 ))|cfa>+^ ( 2 + ^ ) ^ £ w\f(w( Xl +i X2 ))\dw 

(2.7) 

Let us consider the case when x is outside cone C(7ra/2). It is known that, as |x| — > oo 
and x G R 2 \ C(7ro;/2) Mittag-Leffler's function and the partial derivatives have the 
asymptotic form (see [1, 7, 14]): 

+ "»> = - ( , 1+ 1 fa2 )r(i 1 - a) +0 (^) < 2 - 8 ' 

and 

9 .{^(xx + *x 2 )} = -A { L_ }_ j_ + O . (2.9) 



<9xj " 1 dxj (xi + ix 2 ) T(l — a) 



x 



These asymptotics are valid uniformly in the region {x G R 2 \ C(ira/2 + e) | R < \x\} for 
given 7r — na/2 > e > and some Rq » 1. 

Proposition 2.1. Let x G R 2 \ C(na/2). We have, as r — > oo 

fpk( \ (k\x\Y Ji(k\x\) 1 logr /1\ 

^ T) = i^"J x 1+ ix 2 r(l-a)— + ° It J (2 ' 10) 

and 

^ ;T) » {(M\ 2 Mmi) i !2ii + fiv (2 . u) 



<9xj a 9xj I \ 2 y Xi + ix 2 J T(l — a) r \r 



These asymptotics are valid uniformly in the region {x G R 2 \ C(na/2 + e)\R 1 < |x| < 
R} for given n — na/2 > e > and R > 0. 

Proof. From (2.8) and (2.9) we obtain 

- f E a (w( Xl + ix 2 ))dw = -L 1 logr + O (-) ; (2.12) 

r Jo xi + ^x 2 r(l — a) r \r/ 
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jT w\E a {w{ Xl + ix 2 ))\dw = o(^j; (2.13) 
J\ 2 \E' a {w{x 1 + ix 2 ))\dw = o(^j. (2.14) 



J_ r , . i 

r 2 Jo 

Then applying Lemma 2.1 to the case when f(z) = E a (z), from (2.4), (2.5), (2.12) and 
(2.13), we obtain (2.10). Next from (2.6), (2.7), (2.8), (2.9), (2.13) and (2.14), we have 



9 -E k a (x; r) = rV-'E' a {r{^ + *«*)) - f^V — T ^-E a (r(x 1 + ix 2 )) 



dxj ay ' ' ay y " \ 2 / ^1 + ^2 

XjJ 2 (k\x\) \~J o E a (w(xi + ix 2 ))dw + O ^-^ . 



t') 2 | 


(-ir 








- ix 2 


fc|a;|y < 


f (-ir 






[ Xl 


- ix 2 



+m 2 



Then from the equation 

a i (my Mm) 1 _ (*\*\ y f (-■) , - i . ii(tW) _ f *V IjA(tW) l 1 



I \ 2 y + ix 2 I \ 2 / I xi — ix 2 V 2 / [ x± + ix 2 

and (2.12) one obtains (2.11). 
□ 

Next consider the case when x is inside the cone C y (na/2). From (2.2) we have the 
expression 

E k a (x; t) = E a (r( Xl + ix 2 )) - (^pj £ E a (rt( Xl + tx 2 ))J 1 (k\x\VT^t)dt (2.15) 

and a direct computation yields 
d 



dxj 



E k a (x;r) 



= rV~ l ^E' a {r{ Xl +ix 2 )) - (^pj £ E' a (rt(x 1 + ix 2 ))J 1 (k\x\VT^t)dt 

+ri i-i (^M^ £ E' a (rt( Xl + ix 2 )){\ - t) Ji(fc|x|v^7)dt (2.16) 

/ EJrtixx + ix 2 ))Ji(k\x\Vl - t)dt 
Jo 

x\ 2 Xj [ E a (rt(xi + ix 2 ))(l - t)J 2 (k\x\Vl - t)dt. 
Jo 



'kV 
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It is known that there exists a positive constant C such that, for all z G C(na/2) \ {0} 
the estimates 

1 ^V<* C 

\ e «( z ) - - e i < 7, r, 2 M/ 2 ( 2 - 17 ) 



a 



(i + M 2 ) 1/2 



and 

'sP«W-a e } I£ TT^F' (2 ' 18) 

are valid (see [7, 14]). Note that in (2.18) there is no restriction on z in a neighbourhood 
of 0. This is because of < a < 1. 

Proposition 2.2. Given R > and e > let x satisfy R" 1 < \x\ < R and Re(x 1 + 
ix 2 ) l,a > e. Then, as t — > oo we have two formulae: 

^(x;r) = y /< "(- 1 + ^) 1/< "{l + (JL)}; (2.19) 

"^ T > = l:{? T "° ( * I+fai) "°}{ 1+0 (^)}- < 2 -> 



dxj i , 



Proof. Since (2.19) is a direct consequence of (2.15), (2.17) and the estimate 

jf 1 S a (rt(x! + ix 2 ))J 1 (^|x|v / n r I)^ = ^ rl/Q(Xl + ^ 2 ) 1/Q (-^) , (2.21) 



first we give the proof of (2.21). 
Write 

f 1 E a (Tt( Xl + tx 2 ))J l (k\x\VT Zr t)dt = f 1 -er 1,a t 1,a (xi + ^Y /a j l (k\x\VT :r t)dt 
Jo Jo at 

+ £ [E a (rt( Xl + ix 2 )) - ± e r 1/a t 1/a (x 1+ ix 2 y^ j l{k \ x \^T—t)dt 

= 1 + 11. 

Then, from (2.17) we have 

W < C f « 4f« * = „(!il) ( 2 .22) 
Jo (l + r 2 |x| 2 ) 1/2 r\x\ Jo (1 + S 2 )V2 V r / 

provided R^ 1 < \x\ < R. 
Let < 5 < 1. Write 

JO ct 

+ f 1 V^ 7 ^ + - 2 ) Vq j l( Jb|x|v^)* (2 ' 23) 
j<5 a 



= /! + /: 



2- 
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Then, one has 

\h\< i e T 1/a 5^Re(x 1 + ix 2 y/^ (2 24) 



a 



Write 



a Js 



Change of a variable 1 — t 1 ^ = s yields 



J <5 



Jo (1 - s) 1 ^ 

and this gives 

| f 1 e~T 1/a (l ~ * lA *)(*i + ix 2 ) 1/a j l{klx \VT-t)dt\ 
Js 

< a [ l - Sl ' a e -r^sRe ( Xl + ix 2 yl° ds 
~ Jo ' (l-s) 1 "" 

< a f 1 ^" e -r^sRe ( Xl + ix 2 ) x ' a ds 

S { l_(l_^ )} i-i 

a 



< 



5 {1 - a)/a T 1/a Re{x 1 +ix 2 ) lla ' 
Therefore we obtain 

e r l/a Re{x l + ix 2 ) 1/a 

|/2 ' " ^-VV^Re^+u^) 1 /"' (2 ' 25) 
From (2.23), (2.24) and (2.25) one concludes that 

Iae-T 1/a (^ + ^) 1/a = (-^) (2.26) 

provided Re (x x + ix 2 ) 1/a > e. A combination of (2.22) and (2.26) yields (2.21). 
Next we prove: for m — 0, 1 

T ^ Q (rt(xi + ix 2 )){\ - t) m J 1 (k\x\VT^t)dt 
Jo 

(2.27) 

= \e rl/a ^ + «2) 1/a T d-«)/« (a . 1 + ix2 )(i-«)/« f ^ 



a" \r 1/a / ' 
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We have 

JO 

= e r 1 / a (x 1 + tx 2 ) 1 ^ T (i-a )/a{xi + ia . a) (i-«)/« 

JO 

Since 

| e-^i 1 - tl ' a )^ + ^ 2 ) 1/a t (i-«)/ Q(1 _ t ) m Uk\x\Vl^t)dt\ 
Jo 

JO JO 

< 2 

- r 1/a Re(xi +ix 2 ) 1/Q ' 

we obtain 

e r 1 ' a t l l a {x l + tx 2 ) 1 / a T (i- a) / at (i- a )/ a{xi + ix2 )(i-«)/a (1 _ t ) m JMA^t)dt 

Jo 

= e rl/a ^ + 1A y (;ri + ^ 2 )d-)/«o (-L) . 

(2.28) 

Now from (2.18) and (2.28) one obtains (2.27). Similarly, for j — 1,2 we have 

f 1 E a (rt( Xl + ix 2 ))(l - t) m J m+1 {k\AVl^t)dt 
Jo 

(2.29) 

= Ti ^}_ e T lla {x, + tX 2 ) 1 /\(l- a )/ a{xi + ix2) d-a)/a ( M 

provided Re (xi + ia^) 1 / > e and < |x| < R. Note that this is a 'rough' estimate. 

Now from (2.16), (2.18), (2.28) and (2.29) we obtain (2.20). 

□ 

As a corollary of Propositions 2.1 and 2.2 we have immediately 

Corollary 2.1. We have: 

for any regular C 2 curve c with c C C(7ro;/2) 

t 1™ > II^(-;t)IcIIl»(c); — ► °° 

/or any non empty bounded open set U 0/R 2 with U C~R 2 \ C(na/2) 

T \imjE k a (.-,r) ||hi (io — 0. 
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3 Construction of the density 

Definition 3.1. Given uj = (cui, uj 2 ) G S 1 set uj 1 - = (— uj 2 , uJi). Define the function 
E 1 / n (x; s, k, uj) by the formula 

E 1/n (x; s, k, uj) = Ey n {{x ■ uj, x ■ u 1 -); , s > 0. 

From (2.2) we have already known that the function Ei/ n (x; s, k, uj) of x G R 2 satisfies 
the Helmholtz equation Av + k 2 v = in R 2 . Since x ■ {uj + iuj ± ) = {uj\ — iuj 2 ){x\ + ix 2 ), 
the function E\/ n {x; s, k, uj) coincides with the Vekua transform of the harmonic function 
Ei/ n {s{uoi - iuj 2 ){xi + ix 2 )/2)) in R 2 . 

The aim of this section is to construct a density g G L 2 (S' 1 ) explicitly such that 

Hg{x) E 1/n {x; s,k,uj), x G B 2R . 
The starting point is the following fact. 

Proposition 3.1([15]). The Vekua transform of the harmonic function 

^kTp{x x + ix 2 )/2 _|_ e ik(p{x l - ix 2 )/2 _ ^ 

coincides with e lfcr '^. 

Let T be a non empty open subset of S 1 . Given g G L 2 (S' 1 ) the function 

^ e ikJp{x 1 + ix 2 )/2 + e ik<p( Xl - ix 2 )/2 _ l j g ^) d(y ^) 

is harmonic in the whole plane. As a corollary of Proposition 3.1 one knows that the 
Vekua transform of this harmonic function coincides with the Herglotz wave function Hg. 
Taking account the fact mentioned above and the definition of Ei/ n (x; s, k,uj), it suffices 
to construct g in such a way that 

j e i*p(*i + «2)/2 + e iM*i - <x 2 )/2 _ ! j „ El/n ^ZJ{ Xl + «r 2 )) (3.1) 

where uj = uj 1 — iuj 2 . 

Using the power series expansion of Mittag-Leffier's function, one knows that if g satisfies 
the system of equations 

r(JTT)(l)7 r «^' = ' m = 1 -- < 3 ' 2 > 

and 

1 / ik\ m r 1 /s\ m 

' 1 j^ gi(p)da{(p)= (_-^ ?m = ,l,..., (3.3) 



T(m+ 1) v 2 



r(- + i 

n 



then (7 satisfies (3.1) exactly. Now consider the case when T = S 1 . We construct g in the 
form 

oo oo 

g{tp) = £ + £ /?_ m p™. 

m=0 m=l 



13 



Since 



f <p m g(cp)da{(p) =p m , f (p m g(cp)da{(p) = /?_ m , 



/5 1 

from (3.2) and (3.3) we get /3_ m = 0, m — 1, 2, • • • and 

1 r(m + l) /scJ\ m 



V n 



Then g becomes 

» 1 T(m+1) /sZJy9\ m / „x 

^ ) = £w^^bd (3 - 4) 



m=0 r ( — + 1 

n 



This is always divergent. So we consider a truncation of (3.4): 

?(m + 1) / suJc 

(n +1 ) 



1 ^ r(m + l) fs57^ 
ffJV (^; s, fc, w) = — ^ — ^ v ( — ) (3.5) 



m=0 r 

where N — 1, • • •. Then one obtains 

^ | e iM^i + ^ 2 )/2 + e ik<p( Xl - tx 2 )/2 _ ^ s> fc> 



m>nN r 



(3.6) 



This shows g^( ■ ; s, k,u) satisfies (3.1) in this sense. Taking the Vekua transform of the 
both sides of (3.6) we obtain the equation 

Hg N (-;s,k,u)(x)-E 1/n (x;s,k,u) = - £ L \ t^) -U^e*"* (3.7) 

m >nN T hi' 



V n 

where x = (rcos 9, rsin 0). Note that this can be checked also directly and the equation 
Hg N ( ■ ; s,k, uj) (0) = £i/ n (0; s, fc, a;) = 1 holds. 

For our purpose we have to consider how to choose s depending on N. One answer to 
this question is the following and it is the main result of this section. 

Theorem 3.1. Let 7 be the unique positive solution of the equation logi + t/e = 0. 

Let 7 satisfy < 7 < 70. Let {s(N)}n=i,--- be an arbitrary sequence of positive numbers 
satisfying, as N — > 00 



(Rs(N)) n = In + 0(1). 
e 



Then we have, as N — ► 00 



sup \Hg N ( ■ ; s(N), k,u)(x) - E 1/n (x; s(N), k,u)\ 

\x\<2R 

+ sup \V{Hg N (-;s(N),k,u)(x) - E 1/n (x;s(N),k,u)}\ 

\x\<2R (3.8) 

o|«^e W( « +l0g7), ]=o(jr-). 
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Proof. Set 



r(m+l) fs _\ m 7 n ^; m6 ™^ ^ 



MM = £ \ m > (-ZJ J m {kr)e™\ E$(z) = £ - ? - 

m >„Arr hi V 7 m=Or hi 



n J \n 
Then from (2.1) we have, for all x with \x\ < 2R 

\S N (x;s)\< E /J S (ifa) ro = {Ei/ n (z) ~ E$(z)}\ z=Ra . (3.9) 

m>niV T I hi 



Moreover, using the recurrence relation 

Tfl Tfl 

J m +i(£r) = — J m {kr) - J' m {kr), J m _i(£x) = — J m {kr) + J' m {kr) 
kr kr 

one has the formulae 

e ie (|; + i~) J m (A;r)e^ = -A;J m+1 (A;r)e^ + ^; 



Then the formulae 



8 e ie ( 8 1 8 \ e- ie (8 18 
+ ?-^ + —r- l-^-i- 



dxi 2 \8r rd9 2 Idr r 89 



( 8 .1 8 \ i, '" ( 8 .18 



(3.10) 



8x2 2 \8r r 89 J 2 \8r r 89 J 

give the estimate 

\-^-S N (x;s)\ + \A-S N (x;s)\ < sj-{E 1/n {z)-E^ n {z)}\ z=Rs 

+k 2 R{E 1/n (z)-E{%(z)}\ z=Rs . 
From the proof of Proposition 3.2 in [16] one has 

n \ \nN+l . n , 

\E l/n (z) - E${z)\ <E^ J 7Y e ' ' ( 3 - n ) 

and 

J n I |n(AT-l)+Z , n| 

|^-{^i/„W-^(z)}| ^nlzr^M TV 6 ' ' ^ 12 ) 
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Note that these are sharper than (3.5) and (3.6) of Proposition 3.2 in [16]. Consider the 
case when z = Rs(N). Then we get, as N — ► oo 



E 



\nN+l 



\Rez n 



-- 1 r I N + 1 + l - 



= 



(Rs(N)) n ( N+1 ^e (Rs ^y 



V 



r jv + i + 



n 



\ 



) 



n \ \n(N—l)+l 

s(N)n\z\ n - l Y,^r- 



( 



-^Y\N+ l - 



= O 



(Rs(N)) 



n(N+l) e (Rs(N)Y 



FN 



n 



Since 



T[N+- 

n 



O 



(N — 1)! ) ' 



from (3.9), (3.10) and (3.11) to (3.13) we obtain 



\S N (x; s(N))\ + \—S N (x; s(N))\ + \—S N (x; s(N))\ 

OX\ 



O 



'N 2 (Rs{N)) 



dx? 



n(N-l) (Rs(N)Y' 



(N-iy. 



Using the Stirling formula, we have 

(N-iy. 



OlN'^e Ve 



N ( - + log 7 



(3.13) 



(3.14) 



(3.15) 



where {^(iV)}^^... is an arbitrary sequence of positive numbers satisfying, as N — > oo 
£(N) = (j/e)N + 0(1) and < 7 < 7 . Now the conclusion follows from (3.7), (3.14) 
and (3.15). 
□ 

From (3.5) we know that -',s,k)N given by (1.1) has the expression 

dl^if, s, k) N = e~ iky ' v g N (ip; s,k,w), ipe S 1 . 

Then from Definition 3.1, (3.8) and the equation Hg^ n ^( ■ ; s,k) n(x) = Hg N ( ■ ;s,k, uS)(x— 
y), we immediately obtain 

Corollary 3.1. Let {s{N)} N=1 ... be same as in Theorem 3.1. Then for any fixed (y,cu) G 
B R x S* 1 we have, as N — > oo 



sup \Hgff u) ( ■ ; s(N), k) N (x) - E 1/n (x - y; s(N), k,u) 



\x\<R 

+ sup \V{Hgll n J- ; s(N), k) N (x) - E 1/n (x - y; s(N), k,u)}\ = 0(N~°°). 

\x\<R K ' 



(3.16) 
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4 Remarks 



Remark 4.1. It should be pointed out that the density (1.1) satisfies 

(^% / ;. ) (-;^^) i2(51) = i- 

Therefore we obtain the relationship between our indicator function and Kirsch's one: 

K(x) < \l£ u) (s(N)) N \. 

This together with Theorem 1.1 explains why K(x) = in a case that the point x 
can be connected with infinity by a straight line without intersecting D. Note that, if 
x G D, then < K(x) < oo, however, from Theorem 1.1 we know that, for all u G S 1 
limjv^oo |/^")(s(JV))iv| = oo. 

Remark 4.2. Theorem 1.1 does not cover the 'critical' case when both C y (u, ir/2n)r\D ^ 
and C y {u, ir/2n)nD = are satisfied. This is coming from a lack of a necessary uniform 
estimate of the function E a (x; s, k, u>). At the present time we do not know what one can 
say about the behaviour of the indicator function as N — > oo in this case. Note that the 
results in [14] and Theorem 1.1 in [16] completely cover this type case. 
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5 Appendix. Proof of Lemma 2.1 

Write 

u(x;r)-f(T(x 1 + ix 2 )) = -(^\ £ + ix 2 ))dt + \ R(x;r) 



where 



R(x;t)= / f(rt(x! + ix 2 )){ Ji(k\x\) - Ji{k\x\^l - t)}dt. 
Jo 

Using the expression 



n 



>-<*i*i^=^S(Sa( 2 J (i - ir ' <ai) 

one obtains 

i t {Uk\x\VT^t)} = U^\ j 2 (k\x\VT^t). 



Then the mean value theorem yields 



1 (k\x\ N 



|J 1 (A;|x| v / r^t)-J 1 (A;|a;|)|<- f-LIj t. 



From this one gets 

\R(x;t)\ < l - (M\ £ \f(rt(x 1 + ix 2 ))\tdt. 
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Now (2.4) and (2.5) are clear. 
Next from (A.l) we have 

2 



JMAVi^t)} = - Xj (i - t)J 2 (k\x\VT^t) 



and this yields 
d 



dxj 



-{u(x; r) - f(r( Xl + ix 2 ))} = - (^j 2 Xj £ f{rt{ Xl + ix 2 ))Ji(k\x\ y/T=T)dt 

^pj rV- x jf* tf'(rt( Xl + ix 2 ))J 1 {k\x\VT^t)dt 

f 1 f(Tt(x! + ix 2 )){\ - t)J 2 (k\x\Vl-t)dt 
Jo 

{Aj{x; r)+R J {x; r)} 



2 /, \ 2 



k\x\_\ k 



k\x 



V 2 

(A.2) 

where 

2x ■ " *■ 
Aj(x;t) = —r-&Ji(k\x\) f(Tt(x 1 + ix 2 ))dt-Ti 3 ~ 1 J 1 (k\x\) / tf'(Tt(x 1 + ix 2 ))dt 

\Ou J " 

Mk\x\) f 1 f(rt( Xl +ix 2 ))dt 
Jo 



+ 2 X > 



and 

2^ ri 
Ixl 2 jo 



Zx ■ f l 

RAx;t) = -—^ / firtixx + ix 2 )){ Ji(fc|x| Vl - t) - Ji{k\x\)}dt 
\x\ Jo 

-rV- 1 J* tf'(rt{x 1 + ix^iMklxlVT^t) - Ji(k\x\)}dt (A3) 
+ (^) xj £ f(rt{ Xl + ix 2 )){{l - tjUklxlVT^t) - J 2 (k\x\)}dt. 



Change of variables and integration by parts yield 

Aj(x; t) = — / t xi + ix 2 )) 

x 1 + %x 2 

I xi — ix 2 \2 ) I r Jo 

Since 

|{J 2 (A;| a ;|v / r^)} = i(^) J 3 (k\x\VT^t), 
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(A4) 



one knows that 

|(i - t)J 2 (k\x\VT=l) - Uk\x\)\ < |i (^fj 2 + 1| t. 

Using this together with (A. 2), (A.3) and (A.4), we obtain (2.6) and (2.7). 
□ 
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